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VO 

The interactions between charged colloidal particles in an electrolyte may be described by usual Debye- 
Hiickel theory provided the source of the electric field is suitably renormalized. For spherical colloids, we 
reconsider and simplify the treatment of the popular proposal put forward by Alexander et al. [J. Chem. Phys. 
80, 5776 (1984)], which has proven efficient in predicting renormalized quantities (charge and salt content). We 
give explicit formulae for the effective charge and describe the most efficient way to apply Alexander's et al. 
renormalization prescription in practice. Particular attention is paid to the definition of the relevant screening 
length, an issue that appears confuse in the literature. 
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^ . I. INTRODUCTION 

The concept of charge renormalization is widely used in theories on effective interactions in charge-stabilized colloidal sus- 
pensions Ql 0] . The basic idea is to consider the highly charged colloidal particle plus parts of the surrounding cloud of 
qq ■ oppositely charged microions as an entity, carrying a renormalized charge Z e g which may be orders of magnitude smaller than 
the actual bare charge Z of the colloidal particle. Replacing Z by Z e g, one can then use simple Debye-Hiickel-like theories to 
t— I , calculate the effective interaction between two such colloidal particles in suspension. The pioneering work on colloidal charge 
renormalization is the paper by Alexander et al. [5] who proposed to calculate effective charges by finding, within the frame- 
work of the Poisson-Boltzmann cell model, the optimal linearized electrostatic potential matching the non-linear one at the cell 
boundary. In |5J, the prescription how to achieve this task, is given essentially in form of a numerical recipe which is often not 
the most convenient way possible. Making use of simple analytical expressions for effective charges recently suggested 
we here describe an efficient way to realize Alexander's et al. prescription for charge renormalization. In the infinite dilution 
limit of a colloid immersed in an infinite sea of electrolyte, explicit and accurate analytical expressions have been obtained only 
, recently |8], when the Debye length is smaller than the colloid size (spherical or cylindrical). 

In order to illustrate the general significance of this celebrated charge renormalization concept, but also to demonstrate how 
broad its potential field of application is, we have selected a few recent studies which are - in one way or the other - concerned 
with Alexander's et al. prescription. Experimentally, a variety of techniques can be used to determine effective charges of 
surfaces and spherical colloidal particles l9l Hoi [Till , electrophoresis measurements being certainly one of the most 

important. Other experiments, focusing on the colloid aggregation behavior lll4Hl5ll . give direct evidence of the important role 
of counterion condensation and charge renormalization in colloidal systems. On the theoretical side, Yukawa-like interaction 
potentials with effective charges from Alexander's et al. prescription are often used as the standard reference curves in calcula- 
tions of effective interactions between charged colloids in suspensions 1H EE USUI, where, most recently, even discrete 
solvent effects 12 lll22l 12311 and density fluctuation effects 1 24] are considered. Charge renormalization and nonlinear screening 
effects are also important when it comes to investigating the phase-behavior of charge-stabilized colloidal suspensions, be it 
the yet unsettled question of a possible gas-liquid phase-coexistence 125l l26l l27l l28l \2 L A 13 Oil . or the solid-liquid phase-behavior 
I3lll32ll . 3D charged colloidal crystals are rheo-optically investigated in 13311 . while defect dynamics and light-induced melting 
in 2D colloidal systems is in the focus of the studies in 1 34] and [35, 36], respectively. Theories dealing with dynamic properties 
of charged colloidal suspensions are presented, for example, in 13711 . Another field of application are complexation problems: 
the complexation of a polyelectrolyte/DNA with spherical macroions ll38l 13^1 . or, vice versa, the complexation of macroions 
with oppositely charged polyelectrolytes H4<X l44l 14211 . Finally, we mention usage of Alexander's et al. work in recent papers on 
microgels 1 43], on the Rayleigh instability of charged droplets 1 44], on charged spherical microemulsion and micellar systems 
£[, and, as a last and more exotic example, in the theory of dusty plasmas H47I1 . We also note that the phenomenological 
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concept of a Stern layer lEill -successfully applied in many studies on adsorption and micellization of ionic surfactants ll49ll - is 
reminiscent of the renormalization procedure under scrutiny here. 

In section|nl the general framework is presented, while Alexander's et al. prescription is revisited in sectionED] Within such 
an approach, a natural screening parameter «pb appears, that plays the role of an inverse Debye length. In the literature, by 
analogy with the behavior of a bulk system, this parameter is often related to the mean (effective) micro-ion density n* in the 
suspension through the familiar expression Kp B = 87rAsn*, where Xb is the Bjerrum length. Such a relation follows from 
electroneutrality in the latter situation, but becomes incorrect in the confined geometry we shall be interested in. Particular 
emphasis will be put on this issue, when the colloidal suspension is dialyzed against a salt reservoir (section llllBL or put in the 
opposite limit of complete de-ionization (section llHO . 



II. CELL-MODEL, POISSON-BOLTZMANN THEORY AND ALEXANDER'S ET AL. PRESCRIPTION 

We consider a fluid of highly charged colloidal spheres having a radius a and carrying each a total charge — Ze (e is the 
elementary charge). These colloids are suspended in a structureless medium of relative (CGS) dielectric constant e and tem- 
perature T (f3 — 1/kT), characterized by the Bjerrum length, Xb — f3e 2 /e . The suspension is dialyzed against a reservoir of 
monovalent salt ions of a given pair concentration c s (the counterions are also assumed monovalent). Due to the fact that the 
colloid compartment is already occupied by the counterions originating from the macroions, the salt ion concentration in the 
reservoir, c s , is always higher than the average salt ion density n s in the colloid compartment, an effect which is known as the 
Donnan-effect 13(1 l5(l 15111 . Note that we not only have n s < c s but also n s /(l — 77) < c s , where 77 is the volume fraction of 
the colloids [the mean density of salt ions in the volume accessible to microions is thus n s /(l — 77) and not n s ]. By contrast, 
the total concentration of microions in the system (i.e. co-ions plus counter-ions) is larger than the total ion-concentration in the 
reservoir. 

The work of Alexander et al. is based on the Poisson-Boltzmann cell model, an approximation which attempts to reduce the 
complicated many particle problem of interacting charged colloids and microions to an effective one-colloid-problem 1 52J] . It 
rests on the observation that at not too low volume fractions the colloids - due to mutual repulsion - arrange their positions such 
that each colloid has a region around it which is void from other colloids and which looks rather similar for different colloids. In 
other words, the Wigner-Seitz (WS) cells around two colloids are comparable in shape and volume. One now assumes that the 
total charge within each cell is exactly zero, that all cells have the same shape, and that one may approximate this shape such 
that it matches the symmetry of the colloid, i.e., spherical cells around spherical colloids. The cell radius R is chosen such that 
7/ = (a/R) 3 equals the volume (or packing) fraction occupied by the colloids. If one neglects interactions between different 
cells, the thermodynamic potential of the whole suspension is equal to the number of cells times the thermodynamic potential of 
one cell. For a recent and more detailed description of the cell-model approximation see Ref. 1 53]. 

In a mean-field Poisson-Boltzmann (PB) approach, the key quantity to calculate is the local electrostatic mean-field potential 
4>(r) (made dimensionless here by multiplication with j3e), which thanks to the cell-model approximation has to be calculated 
in one WS cell only. It is generated by both the fixed colloidal charge density as well as the distributions n±(r) = c s e T ^ r ^ of 
mobile monovalent ions, and follows from the PB equation, which together with the boundary conditions takes the form, 

V 2 (/>(r) = Kj? es sinh</>(r) a<r<R 
ft ■ V0(r) = ZXb/o 2 r = a (1) 

ft ■ V0(r) =0 r = R 

with ft the outward pointing surface normal and the inverse screening length k ics defined in terms of the ionic strength of the 
reservoir: K 2 es = 8ttXbc s - Here we assume at r = a the constant-charge boundary condition, and impose with the second 
boundary condition at r = R the electroneutrality of the cell. As input parameter we have Z, n ICS , a, Xb and R. But inspection 
of eq. ([0 reveals that in fact only three parameters are really independent, which are K les a (reservoir salt concentration), ZX b/cl 
(colloidal charge) and R/a (volume per colloid, i.e., volume fraction). Note that we have also tacitly assumed the electrostatic 
potential to vanish in the reservoir, where n± then becomes equal to c s . The microionic charge density p(r) can be written as 

K 2 

p(r) = n + (r) — n_ (r) = sinh Mr) (2) 

AttXb 



and the microionic particle density as 



pr(r) = n+(r)+n-(r) = -f- cosher) . (3) 



Hence, 

r R 

Z = 4vr / drr 2 p(r) (4) 
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due to the electroneutrality of the cell, while the total number 2N S + Z of microions in the cell is obtained from 

r ii 

2N S + Z = 4tt / dr r 2 p T {r) (5) 

J a 

where N s is the number of pairs of salt ions in the cell. We define the ("nominal") salt concentration in the system as n s — N s /V 
where V — 4ttR 3 /3 is the volume of the WS cell. This volume should not be confused with that accessible to the microions 
Vfrco = 47r(i? 3 — a 3 )/3 = V(l — rf). The "net" salt ion concentration in the system then reads n s /(l — rf), and due to the 
Donnan effect we have n s /(l — 77) < c s . 

Here is the point where we can describe Alexander's et al. recipe for charge renormalization: 

i) solve the full non-linear problem, eq. Q, find the potential at R, and the total microion density, n± = c s e^ R , 

ii) define the inverse screening constant from the microion density at WS boundary, 

Kp B = 4tt\ b (nf + n R j = «4 S cosh <f> R , (6) 



iii) linearize the boundary value problem in eq. Q about <fi = 4>r, determine the potential solution of the linearized PB 
equation such that linear and non-linear solution match up to the second derivative at r = R, 

iv) compute the effective charge Z e g from this solution by integrating the charge density associated with the linear solution 
from a to R. 

Replacing now (Z, K rP . s ) by (Z e g, ftps) in the effective Yukawa pair-potential of the DLVO (Derjaguin-Landau-Verwey- 
Overbeek) theory j3, 



pv(r) = Z$ i \ B [— , (7) 

one has succeeded in retaining an effective pair-interaction of the simple Yukawa form, even in cases where a high value of the 
bare charge Z may violate the condition for the linearization approximation which this interaction potential is based on. While 
the relevance of defining an effective potential as @ may be questioned at not too low density, the effective charge Z e s and 
inverse screening length kpb are defined without ambiguity. It is also noteworthy that as an exact property of the cell model 
under study 1 55], «pb is related to the pressure P through 

PP =-t^-4b- (8) 



III. ALEXANDER AND COLLABORATORS' PRESCRIPTION REVISITED 

In the original paper of Alexander et al [5], the procedure just described was introduced as a numerical recipe, and it was not 
pointed out that most of this scheme can actually be performed analytically. We now describe the more direct way to calculate 
Alexander's et al effective charge which needs as input nothing but <pR from the solution of eq. Q. 

Let us consider the linearized version of eq. Q. We linearize the PB equation about <pji and find 

V 2 0(r) = K 2 les [ 7o + 0(r)] cosh0 fl (9) 

where <f> — <f> — (j)R and 70 = tanh We observe that the inverse screening length in eq. (|9jl is given by Kj? es cosh <f)R which 
is just Kp B from eq. 0, and thus just the value required in the Alexander et al. prescription. This observation shows that 
Alexander's prescription, in essence, amounts to a certain linearization approximation where the potential value about which to 
linearize, is just the potential at the cell boundary (f>R. Linearizing eq. (0 and 0, the microionic charge and particle densities 
now take the form 

P(r) - -^%(7o + 0M) (10) 
= ^%(1+7O0W). (11) 



4 




With eq. l|6), 7o can be rewritten as follows 

7b=tanhfo=yi- . (12) 

Demanding that <j)(r) as well as its derivative are zero at r = R, we make sure that the solutions of the linearized and fully 
non-linear PB equation agree at the cell edge as required in Alexander's et al. scheme. Taken together the boundary value 
problem of the linearized problem now reads, 

V 2 0(r) = 4b(7o + 0W) a<r<R 

£(r) = r = R (13) 

n ■ V0(r) = r = R 

Note that eq. is a two-point boundary value problem, while eq. fl3l represents a one-point boundary value problem. 



A. Effective colloidal charges 



The solution of eq. dl 31 is given by 



0) = 7o 



u- + /-- 



(14) 



with 



f± = hF ^ R - 1 exp(=F^PB-R)- (15) 

Zkpb 

With eq. dl4> we can now easily calculate the effective charge by following eq. (0J and integrating eq. dlOt over the cell volume, 
resulting in the simple formula 



ry 70 

ZeB — 



K PB Aj 



- j(rep B ai? — 1) sinh[ftpB(-R — a)] + kpb{R — a) cosh[«PB(i? — a)]| . (16) 



By comparison with the numerically calculated effective charges from Alexander's et al. original prescription, we have explicitly 
checked that this formula indeed produces Alexander's effective charges. Note that eq. dl6> also follows from Gauss' theorem 
applied at the surface of the colloid. 

In summary, this then is the simple procedure to calculate renormalized charges: 

i) solve eq. Q to obtain ipn, 

ii) calculate k pb = K^ es cosh <j>n, 

iii) insert this into eq. d!6l > to obtain Z c s ■ 

We emphasize that once eq. Q has been solved numerically, no further numerical fitting procedure is required to match the 
electrostatic potential from the linear and the non-linear solutions at r = Rl In addition, the numerical solution of eq. Q is 
much simplified by taking account of the technical points described in appendix|X](see also appendix|E]for the salt free case). 



B. Effective salt concentration 

Alexander's original paper contains also a prescription how to calculate the effective salt concentration n c s s . Again, we here 
can derive a simple formula. Integrating the particle density in eq. dl It over the cell volume accessible to the microions, as in 
eq. (0, we obtain 

2Nf + Z eS = J^-V bee (l - 7o 2 ) + Z off7o (17) 

with Z c h from eq. d!6l >. so that 

2nf = -/%-(!- 7o 2 )(l - V) - Z cS n c (l - 7 o) , (18) 
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FIG. 1 : The inverse screening length k as a function of Alexander's et al. effective charge Z e g [eq. 1161 1 for various reservoir salt concentra- 
tions c s and two different colloid volume fractions [77 — 0.003 in a), and r\ — 0.03 in b)]. Solid line: Kp B from eq. ^6), dashed line: n 2 from 
eq. <2H . dashed-dotted line: Kp B | 70 _o from eq. 1 1 91 . All curves terminate when the effective charge has reached its saturation value. 



where n c — 1 /V is the colloid density in the suspension and = /V. This formula is again successfully checked against 
the numerical values in Alexander's et al. work. If 70 — > 0, eq. Jl 8I > can be rearranged to give 



l PB 



= 4nX B (Z eS n c + 2nf ) 

To -»0 1 - 77 



(19) 



Many studies, applying Alexander's et al. renormalization concept, can be found in literature where kpb is computed not 
from eq. (|6}, as it should, but from 
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4irX B (Z eS n c + 2n s ) 



(20) 



n 2 = 4:Tr\ B (Z cS n c + 2n s ) , (21) 

with Z c ff from Alexander's prescription and with n s being the actual (as opposed to "effective") salt concentration in the system. 
We here learn that this is certainly not the same as the kpb from eq. (|6jl! Indeed, Eqs. j20t and \2\\ rely on two further 
assumptions, namely that 70 can be set to zero and that n s « nf. Fig. Q checks how good these approximations are. We 
restrict the following discussion to low values of the volume fraction 77 so that expressions J20i and \2\\ coincide (the rather 
formal case of larger 77 will be addressed in section UlI C> . For two volume fractions and various reservoir salt concentrations, 
we calculated in Fig. Q the screening factor as a function of Alexander's et al. effective charge obtained from eq. d!6l >. (i) using 
eq. (solid line), (ii) using the formula 47rAs(Z c ffn c + 2n c ®) from eq. dl9l with ?i° ff from eq. Jl 8I > (dashed-dotted line) and 
(iii) using 4irAB(Z c gn c + 2n s ) from eq. (12 1 i with n s from the non-linear calculation, i.e. from eq. Q (dashed line). As Z e s 
goes to zero, n s — > n° ff , n s — > c s and <j) B —* 0, so all three expressions must lead to the same k 2 = 8it\bc s = n 2 cs . To 
be specific, we took typical values of aqueous colloidal systems for a and X B (a = 60 nm, \ B = 0.713 nm). If however one 
wants to be a more general, one has to specify just K les a, ZX B /a and 77 as these are the independent parameters of the problem. 
Multiplying the values of the x-axis with X B /a and those of the y-axis with 4ir\ B /a where \ B /a — 0.713/60 = 0.012, one 
can transform Fig. Q into a plot that is valid for systems characterized by other values of a and X B . The two colloid volume 
fractions considered in Fig. Q are both easily experimentally realizable. Note that the curves terminate at different Z^f which 
is due to the fact that the saturation value of the effective charge depends on the salt content of the suspension |7|]. 

It is evident from the figure that for low volume fraction (77 = 0.003), it is of no consequence if the screening factor is 
calculated from eq. (|6jl or eq. ( 12 ll . and the error is certainly negligible. At higher volume fraction, however, it is seen that 
both approximations involved in taking eq. (12 It instead of eq. - namely, n s « n c s s , and 70 ~ - take effect: both 
formulae, eq. ( 1191 and (12 ll . fail to give the correct value for kpb at low salt (2c s a 3 = 4) and high Z e g, but the agreement 
between eq. d!9> (dashed-dotted line) and kpb (solid line) improves if c s a 3 increases, while eq. ( 12 1 i still remains a rather poor 
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approximation of kpb- This means that at high salt concentration 70 « produces only a small error, while n s w n e s s is 
always a bad approximation at high volume fraction, regardless the value of c s a 3 . We will refer to expression (12 1 i as a "naive" 
inverse screening length. While such an estimation is inappropriate in the context of Alexander and collaborators' scheme, it 
is noteworthy that it naturally arises from a statistical mechanics treatment of electrostatic interactions in colloidal suspensions 
1 58, 59]. This treatment is however performed within a linear theory formalism and therefore discards the non linear effects we 
are interested in the present article. 

Often, n s is known from the experiment, while c s (and thus 70) is not, so that it seems to be difficult to calculate n c s s from 
eq. i ll 8t . However, it is not. One can obtain n s as a function of c s , from the solution of eq. Q by means of eq. Q and 0. In 
cases where only n s is known from the experiment, one can then use this curve to find the c s corresponding to the known n s . In 
fact, it is immaterial whether or not the experiment was actually performed with the system coupled to a particle reservoir. This 
becomes clear from the following consideration. Take a system coupled through a semi-permeable membrane to a salt reservoir 
with salt concentration c s , and allow for some time till the Donnan equilibrium is reached. Then, the salt concentration in the 
system is n s /(l — 77) < c s . Now, replace the semi -permeable wall by a unpenetrable wall and decouple the reservoir. The 
microion-distribution between the colloids and thus the screening factors will not change. In other words, in cases where the 
system is not coupled to a reservoir one can find a reservoir with an appropriately chosen c s which when coupled to the system 
would leave the microion density distribution in the system unaltered. This means that all our considerations presented are also 
valid for experiments in which the system is not coupled to a reservoir. This, of course, is strictly true only at a mean-field level 
of description. Practically, one then has to proceed as follows: (i) start from a trial value for c s and thus for k 2 . cs = 8it\bc s , (ii) 
solve eq. Q, (iii) calculate n s from eq. l|3} and (|5jl, (iv) vary c s and repeat (i) to (iii) until ac s = cf p is found which leads to a 



that equals the n e s xp from the experiment. The pair (n e s 



~ cp ) can now be used in all the formulae presented above. 



We close this section with a rather general remark. There is actually no need to linearize about the potential at the WS cell 
edge as done in eq. d 1 3t . An alternative is suggested by the following observation. Insert eq. l|3} into eq. (|5jl and linearize about 
a potential </>, 



2JV, 



H- 



L-2 



B 



dr r 2 cosh</>(r) 



dr r 2 



cosh <b + sinh 



(22) 



If one now chooses 



An 

Vfrcc 



dr r 2 (f>(r) 



(23) 



then 



2N° 



Z e S = 



AnX, 



-Vfree COsh ( 



(24) 



so that with n 2 = n 2 cs cosh <j> as in eq. (|6} one obtains 



1-77 



4ir\ B (Z eS n c + 2nf ) 



(25) 



In words: linearizing not about the cell edge value of the potential, but about the average value of the potential in the cell 
(eq. d23t '). leads to effective salt concentrations and effective charges which - when used to calculate the effective ionic strength 
- can be directly related to the inverse screeni ng leng th in the way given by eq. d25i . familiar from the Debye-Hiickel theory. 
This linearization scheme is worked out in l28l 130.15611 . but does not correspond to the original proposal of Alexander et al. 



C. Situation without added salt 



In the limit where the system is in osmotic equilibrium with a reservoir of vanishing salt density (c s — > 0), or in the canonical 
situation where no salt is added to the solution, the PB equation takes the form 

V 2 0(r) = -n 2 e-^ r \ (26) 

where // is a prefactor whose value is determined through the electroneutrality constraint. We have 

4b = 47TA B n5 = [i 2 e -< K (27) 




FIG. 2: Ratio of the "naive" inverse screening length k as given in eq. 1281 over kpb defined from the microion density at the WS cell 
boundary, eq. <27t . as a function of Alexander's et al. rescaled effective charge Z c s\b/o, [eq. 1161 1. for various colloid packing fractions rj in 
the salt free case. The curves terminate at different Z^Xs/a since the saturation value of the effective charge depends on the packing fraction 
& 



Equation dl6l . which provides the connection between the relevant screening parameter kpb and the effective charge Z e g, is 
still correct if 70 = 1 Q|. It is then tempting to compare kpb to the counterpart of expression dl9l : 

9 1 

k 2 = ATTX B Z eS n c . (28) 

1 — 7/ 

This comparison is shown in Fig. |2] Over a wide range of packing fractions, the ratio k/kpb deviates from unity, but less than 
30%. At relatively large bare charges, Z and Z c g differ significantly, so that neglect of charge renormalization in computing the 
screening parameter k [substitution of Z c s by Z in eq. J28H badly fails compared to the "exact" kpb- 

It is instructive to investigate analytically a few limiting cases. We have computed kpbR from the solution of the non-linear 
PB. At fixed volume fraction, this quantity is bounded from above by its saturation value obtained when Z — > 00. In the limit 
where the packing fraction rj vanishes, this saturation value np^R is observed to vanish (although very slowly, as — 77 1 / 3 log 77). 
Hence, Kp^a = rj^'^Kp^R also vanishes, and this piece of information allows to linearize the exact relation ( H6L which then 
takes a simple form: 

Zcff— ~ UkpbR) 2 V- 1/3 , (29) 
a 6 

where the symbol ~ stands for "asymptotically equivalent". In other words, we have 

k pb V ~° 47rA_BZ Gff n c , (30) 

so that k from eq. i28i and kpb coincide in this limit. It may be observed in Fig. [2]that r\ = 10 4 is not small enough to get 
the limiting behavior k/kpb = 1, irrespective of Z (and thus Z c s). In the opposite (and academic) limit where 7/ — > 1, (i.e. 
R —> a), we obtain from (1161 

-5- 1. (31) 

This last result somehow illustrates the relevance of the factor 1 — 77 in (1281 . If the microions densities would have been 
defined with respect to the total volume of the cell and not the sub-volume accessible to microions, we would have obtained the 
alternative naive expression for the screening parameter 



(V) 2 = ATr\ B Z eS n c , 



(32) 
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so that 



Alternatively, in the limit of low bare charge Z we get 



— \/Hv 2/3 -V)^0- (33) 



Z cS ^ Z ^ I (l-l)( KpB ar. (34) 



a 3 \T) 

As a consequence, 

— Z =P 1. (35) 

This feature may be observed in Fig[2]and is compatible with the result of eq. ( 13 0> . 

With added salt, similar argument may be put forward to provide an analytical relation between kpb / k 1cs and the density 
for e.g. large dilutions. Such expressions are nevertheless physically less transparent than those derived here for de-ionized 
suspensions, and have been omitted. 



IV. CONCLUSION 



For colloidal spheres, we have reconsidered the original charge renormalization prescription proposed by Alexander et al. 
|01 . The computation of renormalized charge and salt content has been simplified in two respects: a) by the derivation of 
analytical expressions giving effective quantities as a function of parameters that are directly obtained from the solution of the 
non-linear PB problem; b) by converting the initial two-point boundary value non-linear PB problem into a computationally 
more convenient one-point boundary value problem (see appendices lAl and 151. While we have restricted here to spherically 
symmetric polyions, similar considerations may be applied to the cylindrical geometry, relevant, e.g., to understand properties 
of solutions of charged polymers as for example the DNA molecule. 

Once the effective charge and salt density are known, a "naive" screening parameter k may be defined as k 2 = 4ttX b (Z e sn c + 
2n° ff )/(l — ?/), where the factor 1 — 7/ accounts for the fact that a fraction 77 of the Wigner-Seitz cell is not accessible to the 
microions. Strictly speaking, this inverse Debye length does not coincide in general with the relevant screening parameter 
kpb, which has to be defined from the microions density at the WS boundary. In all the cases investigated here, the difference 
between k and kpb was less than 40%. Moreover, in the salt free case, k and kpb have been shown analytically to coincide 
for both low and large packing fractions (irrespective of the charge), and also for vanishing bare charges. This implies that the 
"naive"expression (3P = (Z c sn c + 2n° )/(l — rf) yields a reasonable zeroth order equation of state for the suspension. On 
the other hand, neglect of charge renormalization, which amounts to defining k or f3P through the bare charge and salt density, 
appears to provide an extremely poor approximation for both kpb and the pressure. 

Acknowledgments: We would like to thank Yan Levin and Jure Dobnikar for interesting discussions, and Mario Tamashiro 
for a careful reading of the manuscript. 



APPENDIX A: NUMERICAL PROCEDURE WITH ADDED ELECTROLYTE 

In this appendix, we propose a few Mathematica® lines of code to solve PB equation for a charged sphere in a concentric 
spherical Wigner-Seitz cell. To this end, all distances are rescaled with the diameter a of the colloid, and the charge expressed 
in units of a/Xg. It is convenient to recast the initial two point boundary value problem Q into a one point boundary value 
problem by assigning an a priori value <j>test to the rescaled potential at WS boundary <f>R. For the situation where n lcs a = 2.0 
at a volume fraction r\ = 0.1, the following procedure finds the corresponding solution (with the arbitrary choice (fitest — 0.1): 

K res = 2.0; 77 = 0.1; R = r ] - 1 / 3 ; <j> test = 0.1 ; sol = (Al) 
NDSolve[ { 24! [r] +rcf,"[r} == K 2 Ies r Sinh[#-]], <t>[R] == <j>tesu <t> '[R] ==0}, 
4>, {r, 1, R}, WorkingPrecision — > 20]; 

The potential <f> may then be visualized as a function of r/a with the command 



Plot[Evaluate[(/>[r]/.soZ], {r, 1, R} ] 



(A2) 



9 




FIG. 3: Bare and effective charges as a function of the "guess" boundary potential 4>test- Here, K rcs a = 2.0 and the volume fraction is 10%. 
The arrow indicates the saturation value for cfitest. 



and the bare charge corresponding to the specific choice made for (f>t es t is obtained as the result of 

( Evaluate[ (f>'[l]/.sol] ) [[1]]. (A3) 

With the above parameters, we get 2.55, which corresponds to the value of ZAs/a. 

The limit <j>test — > corresponds to the limit Z — > 0. However, the solution of our one point boundary value problem only 
exists for tfitest < 4>sat, and the limit (fitest — ► 4>sat corresponds to Z — » oo which is equivalent to Z e g — » Z^. Consequently, 
starting from low values of cf>test> the solution associated with a targeted Z is easily found by dichotomy, adjusting iteratively 
the values of cf>test- For any value of <fit es t, i^pb follows from eq. (jfji and Z e g is computed invoking eq. d!6l >. For the parameters 
used in the example iAl\ . the relation between Z, Z c s and cutest is illustrated in Figure[3] 



APPENDIX B: NUMERICAL PROCEDURE WITHOUT ADDED ELECTROLYTE 

In the salt-free case, it is also convenient to rephrase the problem under study as a one point boundary value problem. At fixed 
bare charge, a possibility would be to consider the limit K res — > of a system in contact with a salt reservoir. In this (formally 
correct) limit, the corresponding 4>t es t diverges, and it turns out to be more convenient to impose that 4>(R) = 0. This choice is 
such that fi = kpb with the notations of eq. (I27> . Poisson's equation (126b is then solved once a trial value has been chosen for 
[i. For instance, with r\ = 0.1 and fx = 1.1, the Mathematica® line 

fj, = 1.1; ?7 = 0.1; i? = ?T 1/3 ; CB1) 
sol = NDSolve[{2u'[r] +ru"[r] == fi 2 r Exp [u[r]], u [R] == 0,u'[R] == 0}, 
u, {r, 1, R}, WorkingPrecision — > 20 ]; <p[r.] := (Evaluate[u[r]/.soZ])[[l]]; 

allows to find the reduced potential <\>. The associated bare charge Z\b/cl is computed a posteriori from Gauss' theorem 

-(Evaluate^!.]]) (B2) 

With the example described in eq. JB1 1 . we get 5.217. Here, /i plays a similar role in the resolution as (fitest in appendix lAl 
The solution of eq. (I26> only exists for /i < /^ sat , and when /i — > /^ sat (i.e. when kpb — * Kpe), the bare charge diverges. 
Implementing a PB-like cell problem along the lines described here constitutes a substantial simplification with respect to the 
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traditional route, usually involving a Fortran or C code with a numerical fitting procedure of the linear and non-linear solutions 
of PB theory. 
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